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SUMMARY 


The  Beattie -Bridgeman  equation  of  state  was  used  to  develop  the 
equations  of  several  of  the  thermodynamic  properties  and  flow  process 
correction  factors  for  gases.  The  increase  in  the  specific  heats  due 
to  the  vibration  of  diatomic  molecules  was  included  by  assuming  the 
molecules  to  be  perfect  harmonic  oscillators.  Thermodynamic  and 
flow  process  equations  are  theoretically  developed.  The  particular 
flow  processes  investigated  were  isentropic  expansions  and  flow  through 
normal  shock  waves. 

The  calculation  procedures  for  isentropic  expansions  and  flow 
through  normal  shock  waves  are  included  in  the  appendix.  Also  in¬ 
cluded  are  the  results  of  several  calculations  in  which  air  was  used 
as  the  media  of  flow.  Graphs  of  the  thermodynamic  properties  and 
Beattie -Bridgeman  correction  factors  for  air  are  included  to  provide 
a  general  picture  of  the  effect  of  temperature  and  pressure. 


NOMENCLATURE 


a,  b 9  C)  A.Q,  Bp 
a 

C 

E 

F 

Si’  Si’  Si’  Si>  Ss  ’  St 
G 
h 
1 
! 
k 
K 
L 
M 
P 


Be attie  - Br id geman  constants  defined  in  appendix 
Speed  of  sound,  ft/sec 
Specific  heat,  Btu/lb 
Defined  by  Eq.  (3) 

Defined  by  Eq.  (44) 

Defined  by  Eq.  (38) 

Defined  by  Eq.  (6) 

Defined  by  Eq.  (24) 

Enthalpy,  Btu/lb 
Defined  by  Eq.  (51) 

Defined  by  Eq.  (47) 

20y/(oy~1'> 

Defined  by  Eq.  (45) 

Defined  by  Eq.  (61) 

Mach  number 
Pressure,  psia 
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q 

R 

S 

T 

u 

v 

V 

w 

X 

y 

o 

p 

PRESUBSCRIPTS 


SUBSCRIPTS 

c 

p 

r 

t 

T 

v 

tr 

vib 

DO 


Heat  energy,  Btu/lb 

Gas  Constant,  (ft 3  lfy)  /  (in.  2  °R  lb-mol) 

Entropy,  Btu/lb  °R 

Temperature,  °R 

Internal  energy,  Btu/lb 

Voliime,  ft3 

Velocity,  ft/sec 

Mass  flow,  lb /ft3  sec 

Defined  by  Eq.  (58) 

Specific  heat  ratio 

Characteristic  atomic  vibration  temperature 
in  lowest  energy  state,  °R 

Density,  lb/ft^ 


Perfect  gas  property 


Compression  component 

Constant  pressure 

Rotational  component 

Total  (isentropic  stagnation  condition) 

Constant  temperature 

Constant  volume 

Translational  component 

Vibrational  component 

Infinite  volume  property  (zero  pressure) 


6 


AEDC-TR-57-10 


INTRODUCTION 


For  wide  ranges  of  pressures  and  temperatures,  the  perfect  gas 
equation  of  state  does  not  define  accurately  the  relation  between  specific 
volume,  pressure,  and  temperature  of  a  gas.  Furthermore,  the  ratio 
of  specific  heats  changes  considerably  from  the  perfect  gas  values  when 
the  temperature  and  pressure  of  the  gas  differ  widely  from  standard 
conditions.  However,  as  long  as  the  temperature  remains  between  300° 
and  700°R  and  the  pressure  remains  between  0  and  30  psia  during  a 
thermodynamic  process,  the  perfect  gas  relations  provide  relatively 
accurate  values  of  the  thermodynamic  properties  of  gases. 

The  Gas  Dynamics  Facility,  Arnold  Engineering  Development  Center 
(GDF  - AEDC)  has  wind  tunnel  circuits  capable  of  producing  pressure 
variations  from  near  0  to  2500  psia  and  temperature  variations  from 
about  70°  to  about  2000°R.  Since  test  results  based  on  the  perfect  gas 
relations  and  constant  specific  heats  are  apt  to  be  in  error,  a  method  of 
correction  is  needed. 

The  Beattie -Bridgeman  equation  of  state  and  the  values  of  the  con¬ 
stants  for  several  gases  were  published  in  1928  (see  Ref.  1).  The 
purpose  of  the  equation  was  to  provide  a  means  of  calculating  more 
accurate  values  of  the  thermodynamic  properties  of  actual  gases. 
However,  the  thermodynamic  relations  developed  from  this  equation 
are  extremely  complicated  and  cumbersome  when  compared  with  the 
relations  developed  from  the  perfect-gas  equation  of  state.  The  addi¬ 
tional  work  involved  is  not  ordinarily  justified  by  the  increased  accuracy 
of  the  calculations;  therefore,  the  equation  has  never  been  used  exten¬ 
sively.  J 

The  need  for  greater  accuracy  in  research  work  has  caused  a  re¬ 
newed  effort  to  develop  methods  of  calculating  gas  properties  with  an 
accuracy  in  keeping  with  the  precision  of  experimental  instrumentation. 
This  work  is  submitted  as  a  step  toward  this  end,  and  with  an  imme¬ 
diate  objective  of  permitting  a  more  accurate  calibration  of  the  GDF 
wind  tunnels . 

This  report  is  an  elaboration  of  the  work  reported  in  Ref.  2  and 
makes  several  changes  such  as:  (1)  the  use  of  engineering  units,  (2) 
the  addition  of  two  more  terms  to  the  infinite  series  expressing  density 
ratio,  (3)  the  use  of  more  recently  determined  values  of  the  character¬ 
istic  temperature  of  vibration  for  molecular  nitrogen  and  oxygen,  (4) 
the  use  of  standard  gas  mixture  proportions  for  air  as  reported  in  Ref. 

3,  and  (5)  the  inclusion  of  two  specific  methods  of  solving  flows  through 
normal  shock  waves. 
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DEVELOPMENT  OF  THE  THERMODYNAMIC  EQUATIONS 


DENSITY 


The  Beattie -Bridgeman  equation  of  state  for  an  imperfect  gas  is 
(see  Ref.  1) 


RT  /i _ c_\ 

i 

v  +  B0\ 

(l_b\ 

i 

A0  j 

v2  \  vT1 ) 

1 

rJ 

! 

{  V  ) 

Letting 


ei  =  Bo“ 


Ao  _  c 
RT  T 3 


e,  =  -^- -B0b~ 

2  RT  °  T* 


Bnbc 


(1) 


Expanding  Eq.  (1 )  and  collecting  terms  we  can  obtain 

_  RT  .  ( d  Ao  c\RT  ( Ao°  d  RT  (B0bc\  rj  (2) 

p^—  +  \B°~RT~r')~+\~RTr~BobYr~)~ir  +  \TT~)~ 


(3) 


and  using  the  relation 
1 


P  =  T 


we  can  write  Eq.  (2)  in  the  form 


p  =  pRT  (1  +  e1p  +  +  e3p3) 


p  t  V 


pk 


(4) 


Let  p  be  represented  by  an  infinite  series  of  the  form 
p  =  +  &P  +  &P2  +  &P3  +  •  •  • ) 


(5) 


Obtaining  successive  powers  of  Eq.  (4),  substituting  in  Eq.  (5) 
and  solving  for  the  values  of  through  &  gives  the  following  relations : 
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i 


Si  =  - 


RT 


2  e 2  —  e2 

(RT)2 


Ss  : 


60^63  -  e2 )  -  e3 

(RT)3 


&  = 


6elei  +  3e22  +  1  e2  (2e2  -  3  e2 ) 

W7T 


(6) 


gs  = 


7  [e2e3  -  4e1  ( e2 2  +  e^)  +  6et3  (2e2  -  e /)] 
(RT)S 


6[  10  et2(  3e22  +  2ete3)  +  lie,4  (2  e2  -  5  e2)  -  2  e23  ]  +  e3  (4e3  -  73  et  e2) 
g6  =  (RD* 


Substituting  in  Eq.  (5)  the  perfect  gas  relation 

_ p_ 

Op  -  RT 

we  obtain  the  density  ratio 

~p  =  1  +  gtP  +  giP2  +  g3p3  +  ■  •  • 


(7) 


SPECIFIC  HEAT 

The  internal  energy  of  a  gas  consists  of  the  kinetic  energy  of  the 
molecules  and  the  atoms  within  the  molecules,  and  the  potential  energy 
of  the  molecules.  When  energy  is  added  to  a  gas,  each  of  these  internal 
energy  factors  takes  up  its  portion  according  to  the  characteristics  of 
the  gas  involved. 

The  kinetic  energy  of  the  molecules  is  dependent  upon  the  number 
of  significant  degrees  of  freedom  of  translation  and  rotation,  and  the 
temperature  of  the  gas.  Most  gas  molecules  take  up  an  amount  of 
energy  very  nearly  equal  to  1/2  RT  per  mol  of  gas  for  each  significant 
degree  of  freedom.  Translation  takes  place  along  the  three  axes  of 
space  for  all  temperatures,  even  approaching  absolute  zero.  Thus,  the 
translation  of  the  molecules  provides  three  degrees  of  freedom.  Rota¬ 
tion  about  the  three  axes  of  space  begins  at  some  finite  temperature 
above  absolute  zero,  and  the  rotational  degrees  of  freedom  become  fully 
activated  at  a  slightly  higher  temperature.  The  number  of  axes  about 
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which  rotation  is  significant  and  the  temperature  above  which  full  activa¬ 
tion  is  attained  are  dependent  upon  the  particular  gas  involved.  Once  a 
particular  gas  is  selected  and  the  number  of  significant  degrees  of  free¬ 
dom  determined,  the  molecular  kinetic  energy  becomes  a  function  of  the 
temperature  of  the  gas  alone. 

The  atomic  kinetic  energy  is  dependent  upon  the  characteristics  of 
the  individual  atoms,  the  number  of  atoms  involved,  and  the  tempera¬ 
ture  of  the  gas.  This  energy  is  held  by  the  moving  atoms  as  they  oscil¬ 
late  with  respect  to  one  another  within  the  molecule  and  is  called  the 
vibration  energy.  *  As  with  the  molecular  kinetic  energy,  once  the 
particular  gas  is  selected  and  the  vibration  characteristics  determined, 
Ijthe  atomic  kinetic  energy  becomes  a  function  of  the  temperature  of  the 
jgas  alone. 

The  molecular  potential  energy  is  the  energy  required  to  force 
the  molecules  into  close  association  with  one  another.  In  any  real  gas 
this  energy  is  a  function  of  any  two  of  the  variables  of  state,  p,  v,  and 
T.  The  perfect  gas  equation  of  state  and  its  accompanying  thermody¬ 
namic  relations  consider  only  the  molecular  kinetic  energy;  they  neglect 
both  the  atomic  kinetic  energy  and  the  molecular  potential  energy.  The 
Beattie-Bridgeman  equation  of  state  and  its  accompanying  thermody¬ 
namic  relations  are  based  upon  a  consideration  of  the  molecular  kinetic 
and  potential  energy.  If  the  atoms  involved  approach  perfect  harmonic 
oscillators,  the  atomic  kinetic  energy  can  be  easily  accounted  for  in 
either  the  perfect  gas  or  Beattie-Bridgeman  thermodynamic  relations. 

The  rotation  of  monatomic  molecules  does  not  represent  a  signifi¬ 
cant  amount  of  energy  since  its  moment  of  inertia  is  negligible  about 
the  three  axes  of  space.  Since  there  is  only  one  atom  present  in  the 
molecule,  there  can  be  no  atomic  vibration.  Thus,  the  internal  energy 
of  a  monatomic  gas  consists  of  only  the  molecular  kinetic  energy  of 
translation  and  the  molecular  potential  energy.  From  the  preceding 
discussion  we  can  see  that  the  internal  energy  of  a  monatomic  gas  should 
be  about  3/2  RT  per  mol  when  the  molecular  potential  energy  is  negligible. 

The  diatomic  molecule  presents  a  problem  which  is  considerably 
different  from  the  monatomic  molecule.  The  three  degrees  of  freedom 


*  The  electronic,  ionization,  and  dissociation  energies  will  not  be 
discussed  since  they  become  significant  only  at  temperatures  above 

the  range  considered  in  this  work. 
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of  translation  remain  the  same.  However,  considering  the  molecule  to 
have  a  dumb-bell  like  form,  where  the  atoms  are  spaced  a  short  distance 
apart  and  held  together  by  a  rigid  bar,  the  moments  of  inertia  about  the 
three  axes  of  space  are  no  longer  negligible.  Visualizing  one  axis 
coinciding  with  the  bar  connecting  the  two  atoms,  we  see  that  the  moment 
of  inertia  about  this  axis  is  negligible  when  compared  with  the  moments 
of  inertia  about  the  other  two  axes  of  space.  From  this  analysis  we  see 
that  for  diatomic  molecules  two  degrees  of  freedom  of  rotation  must  be 
added  to  the  three  degrees  of  freedom  of  translation.  Thus,  the  internal 
energy  of  a  diatomic  gas  should  be  about  5/2  RT  when  the  molecular 
potential  energy  and  the  atomic  kinetic  energy  are  negligible. 

The  potential  energy  of  a  gas  approaches  zero  as  the  pressure 
approaches  zero  (or  the  volume  becomes  infinite)  regardless  of  the 
gas  temperature.  Thus,  one  can  conceive  of  a  hypothetical  expansion 
of  a  gas  by  a  constant-temperature  process  until  the  volume  becomes 
infinite.  Then  the  only  internal  energy  remaining  in  the  gas  is  the 
kinetic  energy  of  the  molecules  and  the  kinetic  energy  of  the  atoms. 
Although  not  rigorously  additive,  a  sufficient  degree  of  approximation 
can  be  obtained  by  separating  the  kinetic  energy  into  the  sum  of  its 
three  components  as  follows: 

ll  oo  —  ^tr  +  Wr  +  u  vib  (8) 

(ur  and  « vib  are  not  independent;  see  Ref.  4). 

Since  the  internal  energy  of  a  perfect  gas  considers  the  translational 
and  rotational  energy,  we  can  write  Eq.  (8)  as 


Uoo  —  ^  vib 


(9) 


The  specific  heat  of  a  gas  is  given  by  the  the rmodynam ic  relation 


where 


(10) 


itq  =  du  +  pdv 


(ID 


To  determine  the  specific  heat  due  to  vibration,  the  pressure  is  assumed 
zero  and  we  can  write  from  Eqs.  (10)  and  (11) 


C  vib  — 


vib 

dT 


(12) 


Integration  gives 


(13) 
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or 


(14) 


Both  Cp  and  Cv  will  be  increased  by  Cvii,  when  atomic  vibrations  are 
considered. 


Tiie  internal  energy  of  the  gas  at  infinite  volume  is  given  by  Eq.  (8) 
as  «oo.  To  obtain  the  total  internal  energy  one  can  assume  a  hypothetical 
constant  temperature  compression  of  the  gas  to  normal  pressure;  then 
the  total  internal  energy  of  the  gas  is  the  sum  of  the  infinite  volume  in¬ 
ternal  energy  and  the  energy  of  compression.  Thus,  we  can  write 


U  =  l*oo  +  Uc  =  0U  +  It  vib  +  Uc 


(15) 


where  uc  is  a  function  of  any  two  of  the  state  variables  p,  v,  and  T. 
Considering  it  to  be  a  function  of  v  and  T  we  can  write 

i*c  =  u(v,  T)  (16) 


from  which 


(17) 


Since  a  constant  temperature  process  is  to  be  used  for  the  compression, 
Eq.  (17)  can  be  written 


duc 


(18) 


Integrating  Eq.  (18)  we  have 


-  -£(&),  d ■ 

When  the  thermodynamic  identity  (-57)  =  T(jrr)  is  usec*» 
Eq.  (19)  becomes 

fv  pv 

T  (fjr)  dv  -  I  Pdv 


(19) 


(20) 


When  we  expand  Eq.  (2),  differentiate  with  respect  to  T  holding  v  constant, 
and  substitute  into  Eq.  (20),  we  get 


3  Rc  A0a  3  BqRc  3  B0bRc 

~  ~  +  ~7T2 


dv 


12 


( 


AE  DC-TR-57- 10 


Integration  and  the  substitution  of  p  =  -i-  gives 


uc  =  R 


(A  O  3  C  \  2 

:  (3 BgC 

A0a\  3 

(  B0bc  \" 

-p[T  +t)  -p 

\2T2 

Tr)+p 

\T2  j 

(21) 


Using  Eqs„  (14),  (21),  and  the  perfect  gas  relation  oU  =  0CvT 
in  Eq.  (15)  produces  the  equation  for  the  total  internal  energy  of  a 
Beattie -Bridgeman  gas  considering  the  contribution  of  the  atomic 
vibrations.  This  equation  is 


u  =  T 


tT 

r 

- 

T  -< 

oCv  +  I  C  vib  dT  +  R 

n(A°  J.  3C>) 

P  \RT  j.%  ) 

2(  3B0c  A0a  \  3  /  Babc  \ 

P  \2T3  2 RT  )+  P  \  T3  )  \ 

Jo 

(22) 


If  we  use  the  perfect  gas  relation 


R  —  oCv{oy  —  1) 


and  let 

Eq.  (22)  becomes 

u=  qCvt\ l  +  (0y-l) 


G  = 


*£ 


(23) 

C  vib  dT 

(24) 

n  ( A-o  3c  \  2  (3Boc  A0a  \  3(^o^c\ 

G-P  (W  +  jr)  - P  {-Jfr  -  m)  +  P  (—) 


(25) 


Using  Eqs.  (10)  and  (11 )  to  obtain  the  specific  heat  for  a  constant 
volume  process  gives 


_ ( du\ 
J v~\dT ) 


(26) 


Substituting  Eq.  (21 )  into  Eq.  (15)  and  differentiating  with  respect  to  T 
with  v  constant  produces 


r  —  (  vib  6  Rc  ( l 

Lv~  \  dT /  dT  dT  T3  9  V 


+  _p 


~3~  P>) 


(27) 


Using  Eq.  (12)  and  the  perfect  gas  relation 


r*  doll 

°  v  =  IT 


we  have  the  specific  heat  at  constant  volume  for  a  Beattie -Bridgeman 
gas  with  the  contribution  of  the  atomic  vibrations  considered 


Cv  —  gCg  +  Cvib  + 


6  Rc 


(1+  T'"  -  "tM 


(28) 
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The  specific  heat  at  constant  pressure  can  be  determined  from  the 
thermodynamic  identity 

C  -C  -  T  (W  df>l  (29) 

^  {dp/dv)rp 

Obtaining  ( dp/dT )2v  and  {dp/dv)  T  from  Eq.  (2)  we  find 


f[ 

_L 

(l  +¥p) 

+  BoP  -  B06p2) 

1  +  2  e  $ 

1  + 

3  e2p2  +4  e3  p3 

(30) 


With  the  equations  for  Cv  and 
in  the  usual  way  as 


Cp,  the  ratio  of  specific  heats  is  determined 


ENTHALPY 


The  enthalpy  of  a  gas  is  defined  by  the  equation 


h  =  U+  JL 


If  we  rewrite  Eq.  (2)  in  the  form 


Tf  =  oCv  T(ay  -  1) 


(«„-  f>)  ♦  o'  (4r  -  B°b 


and  add  Eq.  (25),  we  get  the  enthalpy  relation  for  a  Beattie -Bridgeman 
gas  with  the  contribution  of  atomic  vibrations  considered: 


^  —  oCvT  -i 


oy  +  ( dy  —  1) 


G  +  p  (  B0-  — 0  -  J£-  )  +  p2  ( _  B0b)  +  p3  (  2B°— ) 
V  RT  T*  )  \  2RT  2  T3  /  V  V  T3  > 


Using  the  two  perfect  gas  relations 

r  oGp 

O^V  ~ 


oy 

ah  =  oGp  T 


(31) 


we  obtain  the  ratio  of  the  Beattie -Bridgeman  enthalpy  to  the  perfect  gas 
enthalpy: 

=  1  +  iez^J± 
w  oh  0y 
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ENTROPY 


Entropy  is  a  state  variable  the  same  as  p,  v,  T,  and  u.  Therefore, 
we  may  obtain  the  change  in  entropy  between  any  two  gas  conditions  by 
proceeding  from  one  condition  to  the  other  by  any  desired  method  or 
path. 

The  thermodynamic  relation  for  entropy  is 

■  T  u  ■  (IH  * 


Integrating  between  point  i  and  o,  we  get 


dT  + 


Holding  T  constant  at  7\  and  evaluating  the  change  in  entropy  involved 
when  v  changes  from  vl  to  vQ  and  designating  this  change  as  (S0-S,)a 
we  have 


(  So  -  Si  )  2  = 


p  &), 

%/  V. 

(r=rt) 


dv 


(33) 


Next,  holding  v  constant  at  va,  evaluating  the  change  in  entropy  involved 
when  T  changes  from  Tt  to  Ta  and  designating  this  change  as  (S0-St)t 
we  have 


(S0-Sl)l  = 


(v  =  vQ) 


dT 


(34) 


Finally,  the  total  change  in  entropy  when  proceeding  from  point  1  to 
point  o  is  given  by  adding  Eqs.  (33)  and  (34): 


S0-Sl=(S0-Sl)t  +(S0-St)} 


i: 


Cd 


dT  + 


f  m, 

J  V. 

1  (  T  --  nr%  \ 


dv 


l(v=v0)  -  '‘(r-r.) 

Replacing  cv  in  Eq.  (34)  by  Eq.  (28),  we  get 
T0  f>Tn  /»  To 


B0  2  B 


ob\ 

3  / 


J%T°  fT°  f  T° 

- :fdT+l  dI+1  P  (l  +  p 

r‘  Jt 1  Jti(p-Po) 

Using  Eq.  (23)  and  evaluating  the  first  integral  of  Eq.  (35)  gives 


l 


(35) 
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Let  the  second  integral  of  Eq.  (35)  be  defined  as 


/’ 

Jt , 


dT  =  R  In  |2- 

r  £, 


(36) 


Evaluating  the  third  integral  of  Eq.  (35)  we  have 


f 

Jt,. 


R  i  F'(1  +  ,4-',x)<r-i [*•  (2+^» -nr)  irr - 7?) 

l(p=Po) 

Replacing  these  three  integrals  in  Eq.  (35)  we  have 


(So-S,),  =  R 


In 


(If)""©  *~  {*  +  '*-*  *?-)(*!-*)}  (37) 


Obtaining  from  Eq.  (2)  and  substituting  into  Eq.  (33)  gives 

»-*>■  -/  s  [f  ♦  7  (a.  ♦  ff)  +  -  *.*)  -  J,  (^)  1  *, 

9/v^(T=Tl)  J 


Integrating,  replacing  by  p  and  adding  to  Eq.  (37)  gives 


S0  -  S,  =  R  |  Zn 


(8»  +  ©  -**  (#-■¥) 


Letting 

Zn  F 


3  2 B0bc 

2 

/  S0C 

Bab\ 

3  2Babc 

+  P°  3  T0* 

■  +  Pi 

+  Pi 

V  r,3 

2  ) 

Pl  3  j\3 

=  p  (bo  .+ 

2c  \ 

r3  / 

.  .a  (  fioC 

p  \  r3 

¥) 

-  p3  ■ 

2B0bc 

3  r3 

(38) 


we  have  the  change  in  entropy  between  condition  o  and  condition  l  for  a 
Beattie -Bridgeman  gas  considering  the  contribution  of  atomic  vibrations: 


Eg  Pi  F, 
Ei  po  E0 


(39) 
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DIATOMIC  VIBRATION/ CORRECTIONS 


Diatomic  molecules  have  only  one  normal  mode  of  vibration.  If  the 
vibration  characteristics  of  the  molecule  approach  the  characteristics 
of  a  perfect  harmonic  oscillator,  the  vibration  component  of  specific  heat 
is  given  as  (see  Ref.  5) 


Cvib  =  R 


d/2  T 


sinh  (  6/  2  T  ) 


'  - «  (40> 


When  CVib  in  Eq.  (14)  is  replaced  by  Eq.  (40),  it  becomes 

bT, 


itvib  —  A 


/’ 


en  t 


sinh {6/2  T) 


dT 


To  integrate  let 


Then 


2  T 


dY  =  -  2T*  dT 


V 

'/ 


Note  that  as  T'-O,  y-* «  and  as  T 7,  y-*y.  Thus, 


u  vib  =  — 


Dividing  through  by  RT  we  have 


u  vib 


£ 


dy 


sinh2  y 


=  G  = 

RT  RT 


L  fTCrtt,  JT  - 
Jo 


which,  when  integrated,  yields 

G  =  _ 


2  y 


(41) 


(42) 


exp  (2y)  —  1 

Equation  (42)  gives  the  contribution  of  the  vibration  of  the  atoms  to 
the  internal  energy  of  a  diatomic  gas  provided  the  vibration  character¬ 
istics  of  the  molecule  approach  the  characteristics  of  a  perfect  harmonic 
oscillator. 


Evaluating  Eq.  (36)  in  a  similar  manner,  we  obtain 

»  T 
1  o 


JmT0  PT0  fT> 

Ti  dT-J0  ^ '''Jo  ^ 


dT  =  R  In  4r- 
Ei 


We  find  the  solution  of  the  first  right  hand  integral  to  be 


£ 


C  vib 


dT  =  -R 


£ 


ydy 


sinh2 


=  R 


y  coth  y  —  In  |  sinh  y 


(43) 
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Since  we  are  dealing  with  only  real  values  of  y,  sinh  y  is  always  positive. 
Thus,  |  sinh  y  |  =  sinh  y  .  Substituting  the  limits  we  have 


R  InE 


y  coth  y  —  In  (2  sinhy) 


Solving  Eq.  (44)  for  E  gives 


(44) 


V"'" 

E  = 


exp  [2y/(exp  2 y  -  l)] 
1  — exp(—  2  y) 


DEVELOPMENT  OF  THE  FLOW  PROCESS  EQUATIONS 


ISENTROPIC  FLOW 


Since  an  isentropic  process  involves  no  change  in  entropy,  we  may 
equate  the  right  side  of  Eq.  (39)  to  zero. 

Thus,  r 


l 

/ t\°y~1  e  Pi  Fi 
\Tj  E,  P  F 


i  =  / 1  yy-1  e/f  i/p 

\Tl  )  El/Fl  1/Pl 

Multiplying  both  sides  by  p/p,  and  rearranging 

oY 

P  =  (  T  \  °y  _  1  E/F  p/pRT 
Pi  \  7Tl  /  Ei/Fi 


Using  Eq.  (4)  and  letting 

/  F 

K  =  -p  (1  +  e,o  +  e2p2  +  e3p3) 

we  obtain 


Let 


6Y_ 

p  _  /  t\cY-i 
Pi  \T,t  K, 


(45) 


(46) 

(47) 
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as  defined  by  Eq.  (32);  then  we  can  obtain 


h  =  oCpT/ 

Using  the  thermodynamic  relation 

V*  =  2  Un-h) 

and  replacing  h  by  Eq.  (48)  we  have 

V.2oCpT  (2±h--j) 

The  thermodynamic  relation  for  the  speed  of  sound  is 

■'  -y(w)T 

Differentiating  Eq.  (4)  and  substituting  in  the  above  equation  we  get 

a2  =  yRT  [  1  +  2 e,p  +  3  e2p2  +  4e3p3] 


(48) 

(49) 

(50) 


Letting 

we  obtain 


I  =  —  [l  +  2e3  p  +  3e2p2  +  4e3p3] 
oY 


<Z 


a 2  =  /  0y/?r  =  7(0a2 

d  usir 


,  =•  iXJ-RT  '  a  -•*  ^ 

Dividing  Eq.  (50)  by  Eq.  (52)  and  using  Eqs.  (23)  and  (31)  gives 


\  -  i ,  *" 
f,tf  1  ■ 


M2  =  V 

a 


( oy  —  l)  / 


Solving  for  produces 


rp  -PX-1  «#*  +  / 

Tt  2 


/« 


or  solving  for  r  we  have 

T  = 


/t 


(51) 

(52) 

(53) 

(54) 

(55) 


oY  ~  1 


IM 2  +  / 


Rewriting  Eq.  (46 )  in  terms  of  stagnation  quantities  and  substituting 
Eq.  (54)  gives 


oY 


El 

p 


IM 2  +/  \°y_1 


Kt_ 

K 


(56) 
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Solving  Eq.  (56 )  for  M  we  obtain 


M 


oY  ~  1 

'•(#)  W 


(57) 


oy- 1 


i 


Substituting  Eq.  (48)  into  Eq.  (49),  dividing  by  Eq.  (52)  and  using 
the  relations  in  Eqs.  (23)  and  (31)  and  M2  =  -^-we  find  that 

a 


2  ht 
V 1 


M  + 


2  / 


oy  - 1  / 


Letting 


X  = 


M2 


2  **  +  -2— 

oy  -  i  i 


and  solving  for  M2  we  get  the  equation 


U2  = 


2  /  X 


oy-  i  /  l  - ■* 

Replacing  M2  in  Eq.  (55)  by  Eq.  (59)  we  obtain  the  solution 

T  =  rt/f  (1-X) 


(58) 

(59) 

(60) 


By  substituting  Eq.  (48)  into  Eq.  (49),  multiplying  both  sides  by 
and  using  Eqs.  (23)  and  (31 )  it  is  found  that 


Pv 


Letting 


and 


0Y-  i  (2  ht-V 


p  V  pRTl  20y 


L  = 


(^) 


k  = 


pRTJ 

2  0y 
oy  -  1 


(61) 


we  obtain  the  result, 


p 

Pv 


=  Jl 

A  \  V  ) 


(62) 
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FLOW  THROUGH  NORMAL  SHOCK  WAVES 

In  order  to  determine  the  relation  between  the  variables  of  state  on 
the  forward  and  aft  side  of  a  normal  shock  wave  it  is  necessary  to  apply 
three  equations  to  the  flow  through  the  shock  wave: 


Conservation  of  Mass 


w  =  PlVt  =  p2v2 

(63) 

Conservation  of  Momentum 

1 

11 

al3 

1 

£|3 

(64) 

Conservation  of  Energy 

=  hti 

(65) 

NOTE:  The  subscripts  1  and  2  apply,  respectively,  to  conditions 
forward  and  aft  of  the  shock  wave. 

Substituting  Eq.  (63)  into  Eq.  (64)  and  using  Eq.  (62)  gives 

'2  h. 


2  h.  y  v 

t2  v  1  v  2  l  T 

T /  I  ~ 


V  *  V, 


-  1  =  k  l  1 


V.  I  \  V 

From  Eqs.  (63),  (58),  and  (65)  it  follows  that 


~k) 


Multiplying  through  by  xt  and  collecting  terms,  we  get 

P  2 


/Pi_y 
\P2  ) 


Xt  (Lj  -  k) 


P 1 
P2 


X,  (k  -  £,,)  +  L, 


-L2  =  0 


(66) 


From  Eqs.  (58),  (65),  and  (63)  we  can  establish  that 

£1  =/M2  (67) 

Xi  \pJ 

Also,  from  Eqs.  (48)  and  (65)  we  find 

\  4  -  4  4  <68> 

Let  Eq.  (53)  be  written  in  the  form 

M2 2  i  (J t2  T 

W  “  Wi  Utl  Tti  -  Tt  T\) 
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By  using  Eqs.  (68),  (60),  and  (67)  it  follows  that 


- 


/  T 
i2i  2 


Write  Eq.  (62)  in  the  form 


Pi 

Pi 


p2L2  (1  -  V/2 Atj) 


Pi^i  (1-  F1V2Ati) 
in  which  we  substitute  Eq.  (58)  to  get 


Pi 

Pi 


P  2^ 2  (  1  ~  j^a) 
PiL.Cl  -  AT,) 


(69) 


(70) 


Substituting  Eq.  (63)  into  Eq.  (64),  dividing  through  by  two  times 
Eq.  (48)  and  using  Eq.  (65),  we  obtain 

,  Pi  Pi2 


Pi 


Pi  Vi 2 


Pi 


P'n'n  2AtJ  2oCPTtlJtl 


2o  cp  rtf  2  A 

Using  Eq.  (58)  and  the  relation 

p  = 


2h 


1 1 


9_ _ P_ 

op  ^  T 


Pi 


_ I _  *L- 

*oCpTuJti  Up  A  RT2 


=  Pi 


- I -  +  (£\ 

°cp rti^ti  Up/ i  RTi 


Using  the  perfect  gas  relation 


R 


o^pioY  ~  l) 
oy 


and  solving  for  Pi,  we  get 


Pi  =  P2 


1 _ 1 

e-T 

(£) 

1  x>Ttihi  +  ( 

2  ’ 

'  oy-i  > 
i  2oy  / 

i - 1 

« 

E-< 

1 _ 1 

N 

E-i 

(5) 

1  *iV*i  + 

a 

(Vj 

I - -1 

e-T 

(71) 


\ 
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CONCLUDING  REMARKS 


Numerous  equations  have  been  derived  in  the  preceding  pages.  For 
the  convenience  of  readers  who  are  interested  in  the  basic  equations 


required  for  the  calculation  of  specific 
is  provided: 

gas  properties, 

the  following  list 

Property 

Eq. 

Property 

Eq. 

CP 

30 

K 

45 

Cv 

28 

L 

61 

p/ Op 

7 

E 

44 

l 

51 

F 

38 

J 

47 

S0  -"Si 

39 
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APPENDIXES 

A.  CONSTANTS  FOR  AIR 

B.  CALCULATION  PROCEDURES 

The  calculation  of  isentropic  expansions  of  air  by  the  method  pre¬ 
sented  is  both  simple  and  rapid.  Usually  the  second  iteration  provides 
a  sufficiently  accurate  result.  The  difference  between  the  results 
obtained  by  using  the  perfect  gas  equations  and  the  Beattie  - Br idgeman 
equations  is  generally  less  than  5  percent  for  the  range  of  pressures 
and  temperatures  investigated. 

The  calculation  of  air  flows  through  normal  shock  waves,  by  the 
method  presented,  is  rather  tedious  and  time  consuming  if  a  desk  calcu¬ 
lator  is  used.  However,  most  electronic  digital  computers  can  handle 
the  problem  easily  if  it  is  necessary  to  have  a  greater  accuracy  than 
provided  by  the  perfect  gas  equations.  The  difference  between  the 
results  obtained  for  air  with  the  perfect  gas  equations  and  the  Beattie  - 
Bridgeman  equations  is  generally  less  than  5  percent  for  the  range  of 
pressures  and  temperatures  investigated. 

It  is  interesting  to  note  that  the  Beattie  -Bridgeman  calculations 
predict  an  increase  in  stagnation  temperature  across  a  strong  normal 
shock  wave.  Figure  8  of  Ref.  7  shows  a  recovery  factor  slightly 
greater  than  one  at  the  stagnation  point  for  the  Mach  3.  24  condition, 
verifying  this  prediction.  The  fact  that  the  recovery  factor  decreased 
at  Mach  4.  18  and  4.  92  can  probably  be  attributed  to  the  low  accompany¬ 
ing  ph .  Note  that  the  table  of  normal  shock  wave  calculations  shows  a 
decrease  in  TtJ  with  decreasing  ph  for  constant  Tti  and  ptl/pt2  • 
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APPENDIX  A 

CONSTANTS  FOR  AIR 


The  Beattie -Bridgeman  constants  used  in  the  numerical  computations 
for  air  are: 


A0  =  4906.5  ..  ft'1t>  2 

lb  mol  in. 


a  =  0.30931  ■  ■ 

lb  mol 


Ba  =  0.73860 


lb  mol 


b  =  -0.17636  ,  *t3' 

lb  mol 


c  =  4.0543  X  104  -^-0R; 

lb  mol 


Other  constants  used  in  the  computations  are: 
0y  =  1.4000 


oCv  =  0.17141-1^- 


R  =  10.729  op  0.068561  ,***" 

R  lb  mol  lb  u  R 

From  the  composition  of  standard  air  given  in  Ref.  3,  it  is  noted 
that  N2  and  02  are  the  only  polyatomic  molecules  of  significant  concen¬ 
tration.  Hence,  we  may  consider  a  composition  of  0.  78088  N2  and 
0.  2095  02  to  determine  the  vibration  components  of  the  various  thermo¬ 
dynamic  quantities.  Due  to  the  form  of  Eq.  (40)  the  specific  heat  of 
vibration  for  oxygen  and  nitrogen  will  be  additive.  Thus,  we  can  write 


^vib  air  —  0*78088  Cvib^y^  +  0.2095  Cvib£j 


Theta  may  be  calculated  from  the  equation 

~  h  co  q  c 


where 


h  =  Plancks  constant  = 

k  ~  Boltzmanns  constant  = 

C  c:  Speed  of  light  » 

co€  =  Normal  mode  of  vibration  +  c, 


6.6234  x  10  ”27  erg  sec 
1.3803  x  10“16  erg/°K 
2.9978  x  10l°  cm/sec 


(72) 
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Unfortunately,  the  value  of  <ue  for  a  given  molecule  has  several  values 
which  are  dependent  upon  the  state  of  the  molecule.  However,  since  the 
temperature  range  covered  in  this  work  is  too  low  to  shift  an  appreciable 
proportion  of  the  molecules  into  the  higher  electronic  levels,  the  normal 
or  unexcited  state  is  used.  Reference  6  gives  these  values  as 

=  2359.6  -i- 


CO 


JeO,  =  1580.4  ~~ 


Thus,  we  calculate 

9n2  =  6109.7  °R 

60  =  4092.1  °R 


The  specific  heat  of  air  due  to  vibration  can  now  be  calculated  from 
Eqs.  (40)  and  (72)  using  the  above  values  of  0.  Thus,  we  obtain 


^vib . 


706.85/r  T 

245.22 IT 

_  sinh  (3054.9/  T)  J 

_  sinh  (2046.1/  T) _ 

(73) 


In  a  similar  manner  G  for  air  is  determined  from  Eq.  (42)  to  be 

Q  .  =  -  4770.9/  T _ 857.29/  T  (74) 

alr  1  -  exp  (  6109.7/f )  1  -  exp  (4092.1 /T) 

Since  laE  is  desired  for  Eq.  (39),  we  must  write 
lnEail  =  0.78088  In  ENj  +  0.2095  In  E0i 

From  this  we  obtain 

^air  =  £,V/-78088  E0^™* 

By  using  Eqs.  (42)  and  (44)  the  final  result  is  found  to  be 

0.2095 

Enir  =  \, - 1  U - 1  (75) 


exp  GNi  I  °-7»088 

exp  Gq 

_1  —  exp  (— 6109.7/71)  _ 

1  —  exp  (--4092.1/T)_ 
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APPENDIX  B 

CALCULATION  PROCEDURES 


ISENTROPIC  EXPANSION 

After  an  isentropic  expansion  from  the  stagnation  conditions  Ttl  and 
Ph  to  a  pressure  p,,  the  Mach  number,  Mlt  and  temperature,  7,  can  be 
obtained  as  follows: 

(1)  Determine  the  perfect  gas  values  of  the  Mach  number,  0Ml  and 
temperature  07’1. 

(2)  Obtain  /, ,  /,  and  K,  for  the  values  of  p,  and07t.  Obtain  Jtl  and 
Kti  for  the  values  of  pti  and  7  . 

(3)  Solve  Eq.  (57)  for  Ml  and  Eq.  (55)  for  7,. 

(Note:  The  last  value  obtained  for  a  quantity  should  be  substituted  in 
subsequent  equations  which  require  that  quantity.  For  example,  the 
value  of  Mt  obtained  with  Eq.  (57)  should  be  used  in  Eq.  (55)  to  calcu¬ 
late  7,). 

(4)  Obtain  new  values  of  /t,  7t,  and  K%,  using  pt  and  the  value  of  Tlt 
obtained  from  solving  Eq.  (55). 

(5)  Solve  Eq.  (57)  for  and  Eq.  (55)  for  7,.  Repeat  this  iteration 
process  until  successive  approximations  produce  the  same 
values  of  Uy  and  7,  to  the  accuracy  desired. 


CALCULATED  VALUES  FOR  AIR 


Eqs. 

Ptl  Psia 

7,  °R 

1 1 

Pi  psia 

7,  °R 

M, 

PG 

2200 

2200 

0. 150 

142 

8.52 

BB 

155 

8.46 

PG 

2000 

2200 

0. 100 

130 

8.93 

BB 

141 

8.86 

PG 

2000 

1500 

0.750 

158 

6.53 

BB 

162 

6.56 

PG 

2000 

1000 

1.50 

128 

5.84 

BB 

127 

5.89 

PG 

1500 

1000 

10.0 

239 

3.99 

BB 

238 

4.02 
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NORMAL  SHOCK  WAVES 


With  Tti,  ptl  and  p{J  Known 


When  a  total  head  tube  is  placed  in  a  supersonic  stream  it  is  generally 
assumed  that  the  pressure  measuring  instrument  involved  reads  the  total 
pressure,  Pt2 ,  behind  a  normal  shock  wave.  The  stagnation  conditions  ptl 
and  Tti  ahead  of  the  shock  are  generally  known.  It  is  often  desirable  to  calcu 
late  the  free-stream  Mach  number  Mx  when  these  three  quantities  are  known. 
The  following  procedure  will  yield  U1  as  well  as  Pl,  Tt»  p2,  T2>  Tt2,  and 

Pi/ Pa  : 

(1)  Assuming  Tti  =  Tti  and  using  the  values  of  ptl,  and  pt2,  determine 
1  tl  and  Jt2 . 

(2)  Solve  Eq.  (6  8)  for  Tfa. 

(3)  Determine  Jt2  for  pt2  and  the  value  of  r  obtained  from  Eq.  (68). 

(4)  Solve  Eq.  (68)  again  for  Tt2 .  Repeat  this  iteration  process  until 
successive  approximations  of  Tt2  produce  the  same  value  to  the 
accuracy  desired. 

(5)  Calculate  the  perfect  gas  values  of  0/W, ,  0TX,  oPi,  0T2,  oP a  and 


(6)  Determine  /1S  Jlt  (P/oP)1P  for  0Tl  and  oPi  ;  J2,  K2,  (p/0p)2  for 

oT2  and  0p2;  and  Kt2  for  T  t2  and  pt2. 

(7)  Solve  the  following  equations  in  order:  Eq.  (58)  for  Xlt  Eq.  (67)  for 
X2P  Eq.  (60)  for  Tlt  Eq.  (60)  for  T2,  Eq.  (46)  for  p2  where  Eq.  (46) 
is  written  in  the  form: 


Pa 


Pt a  K * 


K 


t2 


and  Eq.  (71)  for  p,. 


(NOTE:  The  last  value  obtained  for  a  quantity  should  be  substituted  in  subse¬ 
quent  equations  which  require  that  quantity.  For  example,  when  calculating 
pt  with  Eq.  (7 1 )  the  values  of  Xlt  X2l  Tlt  T2>  and  p2  obtained  in  the  pre¬ 
ceding  calculations  should  be  used  in  Eq.  (71).  ) 


(8)  Determine  Iit  /,,  Klt  L2  for  the  new  values  of  Tx  and  p,  ;  and  J2, 

K2  ,  L2  for  the  new  values  of  T2  and  Pa. 

(9)  Solve  for  pjp2  and  Mlt  respectively,  from  Eqs.  (66 )  and  (57). 

(10)  Repeat  the  above  procedure  for  Xlt  X2 ,  Pl.,  T1S  p2,  and  T2. 

Determine  new  values  of  /,,  Jlt  Kla  Lx>  J2t  F^,  and  L2t  and  solve 
again  for  pJ Pi  and  Mt,  This  iteration  process  should  be  continued 
until  all  calculated  values  repeat  to  the  accuracy  desired. 
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ptl,  and  Mx  Known 

The  static  conditions  before  and  after  a  normal  shock  wave  and  the 
stagnation  conditions  after  the  shock  can  be  determined  when  Mlf  ptlJ 
and  Tti  are  known.  The  gas  properties  obtained  during  the  calculation 
process  are:  plt  Tx,  p2,  TJt  pt2,  oTx,  and  pJp2.  The  calculation  proce¬ 
dure  is  as  follows: 

(1)  Calculate  the  perfect  gas  values  of  0pt,  0Tx,  0p2,  aT2  and  apt2, 

(2)  Determine  lu  fl,  Kt,  for  0p,  and  07;  Jtl,  Kti  for  ptl  and  Tti. 

(3)  Solve  the  following  equations  in  order:  Eq.  (58)  for  Xu  Eq. 

(60)  for  7,  and  Eq.  (46)  for  p,  in  the  form 


qy 


(Note:  The  last  value  obtained  for  a  quantity  should  be  substituted  in 
subsequent  equations  which  require  that  quantity.  For  example,  when 
calculating  Eq.  (46)  for  p,,  the  value  of  7,  obtained  with  Eq.  (60)  should 
be  used  in  Eq.  (46). 

(4)  Determine  /,,  /„  Kx  for  these  new  values  of  p,  and  Tt. 

(5)  Repeat  the  above  calculation  procedure  for  X2,  7„  and  pt .  This 
iteration  process  should  be  continued  until  all  calculated  values 
repeat  to  the  desired  accuracy. 

(6)  Determine  7,  for  p,  and  Tx;  J2,  L2  for  0p2  and  QT2. 

(7)  Solve  Eq.  (66)  for  (px/p2)  ,  Eq.  (67)  for  X2 ,  Eq.  (60)  for  T2,  and 
Eq.  (70)  for  p2 . 

(8)  Determine  J2,  L2  for  these  new  values  of  p2  and  T2. 

(9)  Repeat  the  above  calculation  procedure  for  ( pjp2 )  ,  X2,  T2,  and 
p2 .  This  iteration  process  should  be  continued  until  all  calcu¬ 
lated  values  repeat  to  the  accuracy  desired. 

(10)  Assuming  0Tti  =  Tti  determine  aJt2  and 

(11)  Solve  Eq.  (68)  for  Tt2  and  Eq.  (46)  for  p(2  in  the  form 

pY 

0°y"1 
Tj 

(12)  Determine  Jt2  and  Kt2  for  these  new  values  of  pt2  and  T t2 . 

(13)  Repeat  the  above  calculation  procedure  for  Tt2  and  p(1 .  This 
iteration  process  should  be  continued  until  all  calculated  values 
repeat  to  the  accuracy  desired. 
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CALCULATED  VALUES 


Eq. 

Pt, 

T 

1  h 

Pt2 

Pi 

L 

Ps 

T2 

Pi/ Pi 

K 

PG* 

4000 

1900 

40 

0.519 

147 

35.9 

1842 

0.  181 

1900 

7.  71 

bb** 

0.508 

154 

36.2 

1877 

0.  171 

1925 

7.  78 

PG 

3000 

1900 

30 

0.389 

147 

27.0 

1842 

0.  181 

1900 

7.71 

BB 

0.385 

155 

27. 1 

1871 

0.  172 

1919 

7.74 

PG 

2500 

1900 

25 

0.324 

147 

22.5 

1842 

0.  181 

1900 

7.71 

BB 

0.  323 

156 

22.6 

1867 

0.  172 

1915 

7.72 

PG 

2000 

1900 

20 

0.259 

147 

18.0 

1842 

0.  181 

1900 

7.  71 

BB 

0.260 

156 

18.  1 

1864 

0.  172 

1912 

7.70 

PG 

1500 

1900 

15 

0. 195 

147 

13.5 

1842 

0. 181 

1900 

7.  71 

BB 

0. 195 

156 

13.  6 

1861 

0.  172 

1909 

7.  68 

PG 

1000 

1900 

10 

0.  130 

147 

8.  98 

1842 

0.  181 

1900 

7.71 

BB 

0. 131 

157 

9.04 

1858 

0.  172 

1906 

7.  66 

PG 

500 

1900 

5 

0.0648 

147 

4.49 

1842 

0.  181 

1900 

7.  71 

BB 

0.0659 

157 

4.  52 

1855 

0.  172 

1903 

7.  64 

PG 

2500 

2200 

7 

0.  0523 

101 

6.  31 

2134 

0.  175 

2200 

10.20 

BB 

0.  0526 

110 

6.36 

2167 

0.  163 

2218 

10.  10 

PG 

2500 

I960 

11 

0.  0997 

108 

9.  91 

1901 

0.  176 

1960 

9.24 

BB 

0.  0994 

115 

9.97 

1928 

0.  167 

1976 

9.23 

PG 

2000 

1900 

15 

0.  171 

131 

13.  5 

1844 

0.  179 

1900 

8.22 

BB 

0.  171 

139 

13.6 

1865 

0.  170 

1912 

8.21 

PG 

1500 

1900 

100 

3.  19 

328 

88.  8 

1837 

0.  201 

1900 

4.  90 

BB 

3.20 

348 

89.4 

1855 

0.  192 

1909 

4.87 

PG 

1000 

1700 

50 

1 . 40 

259 

44.  5 

1645 

0. 197 

1700 

5.28 

BB 

1.38 

271 

44.7 

1657 

0.  189 

1705 

5.26 

^Perfect  gas  equation 

^^Beattie-Brid^eman  equation 
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ILLUSTRATIONS 


Several  of  the  thermodynamic  properties  and  Beattie -Bridgeman 
correction  factors  for  air  have  been  calculated  using  the  equations 
derived  in  the  body  of  this  report.  These  quantities  are  contained  in 
tabular  form  in  Ref.  8.  Graphs  of  the  tabulated  values  are  included 
herein.  Quick  approximate  solutions  of  various  thermodynamic  prob¬ 
lems  can  be  obtained  using  the  graphs.  A  Mollier  diagram  for  air  is 
also  included;  the  values  of  enthalpy  and  entropy  were  not  included  in 
the  tables  of  Ref,  8. 
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SPECIFIC  HEAT  AT  CONSTANT  PRESSURE, C 


Fig.  2.  Specific  Heat  at  Constant  Pressure  vs  Pressure 
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Fig.  3.  Specific  Heat  at  Constant  Volume  vs  Pressure 
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